Tunable polaronic conduction in anatase Ti0 2 
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Oxygen vacancies created in anatase Ti02 by UV photons (80 - 130 eV) provide an effective electron- 
doping mechanism and induce a hitherto unobserved dispersive metallic state. Angle resolved photoemission 
(ARPES) reveals that the quasiparticles are large polarons. These results indicate that anatase can be tuned 
from an insulator to a polaron gas to a weakly correlated metal as a function of doping and clarify the nature of 
conductivity in this material. 



The anatase structural phase of titanium dioxide (Ti02) can 
be the key element in novel applications. Whereas extensive 
work has been focused on its famous photocatalytic behav- 
ior (TJ-O, more and more proposed devices, such as memris- 
tors IH, spintronic devices 0, and photovoltaic cells (SHU, 
rely on its less well-known electronic properties. In partic- 
ular, anatase has been recently suggested as a candidate for 
replacing the In-based technology for transparent conducting 
oxides | 9 ] in a wide range of applications from solar cell ele- 
ments, to light-emitting devices, to flat panels, to touch- screen 
controls ifTOl . The crucial quantity for the figure of merit in 
these devices is conductivity, and it is therefore of major in- 
terest to understand and control the electronic properties of 
pristine and doped anatase. 

Stoichiometric anatase is an insulator with a 3.2 eV band 
gap ifTTl but oxygen vacancies, typically present with con- 
centrations in the 10 17 cm -3 range [Q21 US, create a shal- 
low donor level ~10 meV below the conduction band (CB) 
lfT4l . Since large single crystals became available for trans- 
port studies, a better insight has been gained on the influence 
of these donors on the electronic response of anatase. Above 
~60 K, the electrons thermally excited into the CB give rise 
to metallic-like transport. At lower temperatures, the anoma- 
lous increase of resistivity indicates that the charge carriers 
are not bare electrons but polarons fT4l . i.e., electrons coher- 
ently coupled to a lattice distorsion induced by the Coulomb 
interaction. Understanding the properties of such composite 
particles in anatase is important to better engineer the mate- 
rial for targeted applications, where the low electron mobility 
often represents the overall performance bottleneck. We will 
also demonstrate that, from the point of view of fundamen- 
tal physics, anatase represents an excellent model compound 
to study the behavior of the "rare" large polaron quasiparticles 
(QPs), intermediate between localized small polarons and free 
electrons. 

We performed ARPES measurements on TiC>2 single crys- 
tals (Fig. 1(a)) and thin films grown in situ on insulating 
LaA103 and conducting Nb-doped SrTiC>3 substrates. Clean 
(001) surfaces were prepared as described in Suppl. Inf. The 
results presented have been obtained consistently both for sin- 



gle crystals and thin films, and therefore reflect intrinsic prop- 
erties of the anatase phase, independent of the sample prepa- 
ration method. While oxygen defects are always present to 
some extent after the surface preparation, we have found that 
exposure to UV photons induces a much larger amount of va- 
cancies, and provides a substantial electron doping. Hence, 
we could tune the electron density over more than two orders 
of magnitude by varying the beam intensity and the oxygen 
partial pressure during the ARPES experiment, and explore 
samples with carrier densities in the 10 18 — 10 20 cm -3 range. 

Fig. l(c-d) shows two ARPES constant energy (CE) maps 
measured at the Fermi level in two perpendicular sections of 
momentum space, parallel (c) and perpendicular (d) to the 
(001) surface. The two planes completely define the three- 
dimensional (3D) Fermi surface (FS) of an ellipsoidal elec- 
tron pocket, elongated in the k z direction and centered at the 
T point, center of the 3D Brillouin zone (BZ) (Fig. 1(b)). The 
data are in agreement with theoretical calculations for the bot- 
tom of the conduction band of anatase TiC>2 fl5l , and more 
importantly, they establish the existence of conduction elec- 
trons with a well defined FS in electron-doped anatase. Note 
that the closed FS contour proves the 3D nature of the elec- 
tronic dispersion, in contrast with the two-dimensional char- 
acter of the metallic states observed at the surface of SrTiC>3 
(Hun) andKTa0 3 (HI. 

The ARPES intensity map of Fig. 1(e) illustrates the 
energy-momentum dispersion of the conduction states for a 
sample with electron density n e ~ 3.5 x 10 19 cm -3 , as deter- 
mined from the volume of the electron pocket. It consists of a 
shallow QP band with minimum at —40 meV, which crosses 
the Fermi level E F at hp = ±0.12 A -1 . Remarkably, the 
QP band is followed by a satellite with a similar dispersion at 
~100 meV higher binding energy, and by a broad tail. The 
satellite is more clearly resolved in samples with lower carrier 
densities, e.g. in the spectrum of panel (f) measured on a sam- 
ple with n e ~5x 10 18 cm -3 , where also a weaker replica is 
visible ~100 meV below the first satellite. 

The intensity distribution of Fig. 1(e) reveals the clearly 
dispersive nature of the states, but also a substantial renor- 
malization of the spectral function, incompatible with a sim- 
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FIG. 1. (a) An as-grown anatase single crystal, (b) The BZ of anatase. (c,d) CE maps at Ef (T= 20 K, hv — 85 eV) of electron-doped 
anatase (001) in the k x k y (c) and k x k z (d) planes, respectively. The blue lines outline the boundaries of the 3D BZs. (e) E vs. k dispersion 
of the bottom of the conduction band for a sample with n e — 3.5 x 10 19 cm -3 , (f) ARPES intensity measured at k = Jzf for a sample with 
n e c± 5 x 10 18 cm -3 . The solid line is a Frank-Condon line shape. Voigt peaks of width AE = 90 meV (FWHM) are separated by 108 meV, 
while intensities follow a Poisson distribution, (g-i) Cartoon of the polaron formation induced by the photoemission process, showing the solid 
in its ground state (g) and two possible final states (h,i) of ARPES. 



pie scenario of a metal with weakly interacting electrons. 
A parabolic fit yields an effective mass of m* xy = (0.7 ± 
0.05)m e , where m e is the bare electron mass. A compari- 
son with m xy = 0.42m e from a band structure calculation 
fl9l , yields a mass renormalization (m* y /m xy ) ~ 1.7. In- 
deed, at these carrier densities (n e ~ 10 18 — 10 19 cm -3 ) the 
low energy states appear to fit in the intermediate regime of 
the so-called large polarons. The electron-phonon (e-ph) cou- 
pling induced by the ionic anatase lattice causes the QP dress- 
ing and the mass enhancement, but the polaron wavefunctions 
extend over several lattice constants. We stress here that, al- 
though similar claims are valid for the rutile phase of TiC>2 
l20lL the same ARPES experiment repeated on several (001) 
rutile surfaces failed to show any trace of a metallic edge, even 
after a long exposure to the photon beam. 
The well-known Franck-Condon (FC) scenario for an elec- 
tron coupled to a vibrational mode provides a schematic 
but instructive guideline to interpret the data l2ll . In the 
electron-removal spectrum, the main "zero phonon" peak is 
followed by a progression of vibronic satellites, separated by 
the phonon energy Huq, with peak intensities following the 
Poisson distribution l22l . A FC line shape indeed provides a 
good qualitative description of the spectrum of Fig. 1(f) for a 
phonon energy of Hujq = 108 meV, the energy of a longitudi- 
nal optical (LO) E u phonon mode observed by Raman spec- 
troscopy l23l . Within the photoemission literature, the series 
of distinct satellites is hardly observable, and the distinctive 
sign of a polaronic system is the characteristic "peak-dip- 
hump" spectrum I24ti27ll . Fig. l(g-i) illustrates the physical 
picture underlying the polaron state as measured by ARPES. 



After absorbing a photon (g), the solid is left with a photo- 
hole coherently coupled to a phonon cloud, moving through 
the lattice in its ground state (h), or in one of its vibrational 
excited states (i). The case (h) corresponds to the QP band, 
while (i) corresponds to the satellite replica(s). 

In order to get more quantitative information on the effect 
of e-ph coupling in this (unusual) regime, we consider, as ap- 
propriate for a polar material, a Frohlich polaron model. It is 
characterized by the non-local interaction between the elec- 
tron and an LO phonon branch l22l : 
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FIG. 2. (a) Spectra extracted from the intensity map of Fig. 1(e) in 
the range -0.25 < k x < 0.25 A -1 , (b) Spectral function A(k, E) 
for the Frohlich polaron, including the QP term and the leading / = 1 
incoherent contribution, corresponding to the first phonon satellite. 
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ct and bl create, respectively, an electron with wave vector k 
and a phonon with wave vector q. is the band mass of the 
uncoupled electron and the dimensionless constant a defines 
the e-ph coupling strength. For sufficiently small dopings, 
a single-polaron approach is appropriate, since the effective 
polaron-polaron interaction is weak in this limit. 

The ARPES spectrum is proportional, via dipole matrix el- 
ements, to the one-particle spectral function A(k,E). For a 
Frohlich polaron it takes the form: A(k, E) = A^\k, E) + 
A inc 0) $,E). The coherent QP spectrum A { °\k,E) corre- 
sponds to transitions to an excited state where the lattice re- 
mains unperturbed. It follows the renormalized band disper- 
sion E(°\k) = H 2 k 2 /2m* — fi , where fi is the chemical 
potential. The remaining I > contributions are the inco- 
herent parts of the spectrum, and involve I phonon excitations 
in the final state. The leading I = 1 term can be calculated 
analytically. It exhibits a logarithmic singularity at energy 
E^ (k) = E^ (k) — Hujo, i.e., the same dispersion as the QP 
band, and increased broadening (Suppl. Inf.). The I > 1 terms 
yield a fairly homogeneous background in the energy ranges 
—Ihuoo — fx < E < —IfkjjQ. Fig. 2 compares the experimen- 
tal spectra extracted from the intensity map of Fig. 1(e) with 
the calculated spectral function, including the I = and I = 1 
terms, with Hujq = 108 meV. A Gaussian broadening (85 meV 
full width at half maximum - FWHM) was applied to the the- 
ory. It accounts for the coupling to low-energy phonons and 
other scattering mechanisms, and for the instrumental resolu- 
tion (30 meV). The overall agreement between theory and ex- 
periment is very satisfactory, considering that the theory does 
not include the background tail of the I > 1 terms. 

The above calculation is too simple to provide an accurate 
estimation on the strength of the e-ph coupling. We can as- 
sess it more reliably analyzing the intensity distribution of 
the ARPES signal between the two branches. Quite gener- 
ally, spectral weight is transferred from the QP band to the 
incoherent satellites as the e-ph coupling is increased. The 
coherent fraction Z(kp) of the total spectral weight can be in- 
ferred from the experimental spectrum of Fig. 1(f), yielding 
Z{kp) = 0.36. For this value, diagrammatic quantum Monte 
Carlo simulations of the electron addition spectrum deduce an 
e-ph coupling constant a ~ 2 (28), in a regime of interme- 
diate coupling. Using as a starting point the experimentally 
observed mass renormalization factor of 1.7, the same numer- 
ical calculations predict a ~ 2.5, in fairly good agreement 
and in the same coupling regime. 

The size of the e-ph coupling also has an influence on the 
energy broadening of the photoemission signal. This is il- 
lustrated by the spectra in Fig. 3(a), measured between 20 K 
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FIG. 3. (a) Spectra measured between T = 20 K and T — 300 K on 
a sample with n e — 5 x 10 18 cm -3 , showing a progressive broad- 
ening with increasing temperature. The signal extracted from the 
shallow electron pocket was integrated over a momentum range of 
—kF < k x < kF. 0>) Temperature-dependent intrinsic QP line 
width. Colors of the symbols correspond to those of the spectra in 
(a). The line is a Bloch-Griineisen curve. 



and 300 K on a sample with n e ~ 5 x 10 18 cm -3 . The 
spectral line shape exhibits a temperature dependence besides 
the trivial broadening of the Fermi cutoff. The instrinsic 
(Lorentzian) FWHM AEl of the QP peak, extracted from the 
spectra after removal of an experimental (Gaussian) broaden- 
ing, is shown in Fig. 3(b). AEl exhibits an approximately 
linear T dependence above T = 150 K, and saturates below 
T = 50 — 60 K to a value ~ 55 meV, which includes con- 
tributions from impurity scattering and from the finite photo- 
electron lifetime. The data are well described by the Bloch- 
Griineisen curve modeling resistivity in metals |29l . From its 
high temperature limit AEl = 2tt XksT one can extract the 
mass-enhancement parameter A ~ 0.7, which again yields 
(m* y /m xy ) = 1 + A = 1.7. This A accounts for the elec- 
tron interaction with all the phonons, in particular with the 
low-energy acoustic modes, and should not be confused with 
a introduced above, which embodies the coupling with the 
single LO mode at 108 meV . 

We now turn to the doping dependence of the spectra. Dur- 
ing the ARPES measurement, oxygen vacancies are created 
by the photon beam in a thin layer below the surface. If 
the sample is at the same time exposed to a small O2 par- 
tial pressure, a competing re-oxidation process takes place. A 
dynamic equilibrium is reached, which depends on the pho- 
ton flux at the sample. This offers a unique opportunity to 
study the electronic states of anatase at various electron den- 
sities, without changing ex situ the stoichiometry or adding 
extrinsic impurities. In the following experiment, we kept 
the photon flux on the sample constant, and varied the oxy- 
gen pressure between 2 x 10 -10 and 5 x 10 -8 mbar. The 
corresponding electron densities, estimated from the volume 
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FIG. 4. Evolution of the CB states with doping. The ARPES spectra 
{jiv — 85 eV) were measured at T = 20 K along the dashed line 
(k y = 2tt/o) of Fig. 1(c) for samples with n e — 5 x 10 18 cm -3 
(al), n e = 3 x 10 19 cm -3 (bl), n e = 1 x 10 20 cm -3 (cl), and 
vie — 3.5 x 10 20 cm -3 (dl). The same data are also presented as 
E vs. k x image plots in the bottom panels (a2-d2). A doping scan 
with finer steps is presented as a movie attached to the suppl. info. 



of the electron pocket at the bottom of the CB, varied be- 
tween ~ 5 x 10 20 cm -3 and 10 18 cm -3 , respectively. At 
constant oxygen pressure, data could be collected for several 
hours without appreciable changes. 

Fig. 4 shows spectra of the CB for various doping levels. 
Panel (a) for n e = 5 x 10 18 cm -3 shows a very shallow elec- 
tron pocket and well-defined satellites. At T = 20 K, the 
momentum width Ak of the ARPES spectral function at kp 
gives a QP coherence length I = 1/ Ak = 7 A. Therefore, in 
this low-density limit the polaronic QPs cannot move freely 
for more than ~ 2 unit cells. Increasing the carrier density to 
n e = 3 x 10 19 cm" 3 in (b) and to n e = 1 x 10 20 cm" 3 in (c), 
one observes that : i) the QP band dispersion deviates from a 
parabola, and ii) near the bottom of the band, the QP intensity 
is reduced and the satellite spreads into the background. At 
n e = 3.5 x 10 20 cm" 3 the QP band is visible only near the 
Fermi level crossings at ±fcp. For smaller wave vectors, the 
QP intensity is strongly suppressed, with most of the spectral 
weight spread over a broad energy range (the "hump"). 

The dispersion of Fig. 4(d) in the vicinity of E = huo is 
reminescent of the characteristic kink structure encountered 
in several Fermi liquids with more moderate e-ph coupling. 
The evolution of the CB states with doping indeed suggests 
a breakdown of the single-polaron picture, as charges added 



to the CB progressively screen the e-ph interaction. Polarons, 
which are well defined QPs at low density, eventually lose 
coherence and dissociate into an electron liquid coupled to 
the phonon(s). The spectral weight distribution for high dop- 
ings can be reproduced by standard perturbation theory l22l 
(Suppl. Inf.). 

Our observations shed light on the conduction mechanisms 
taking place in anatase-based devices, and in particular on the 
role of the e-ph coupling, which has been shown to repre- 
sent the dominant scattering process at typical operating tem- 
peratures in pristine and doped films |9 |. The tunability of 
the doping level by UV (or e-beam) illumination over a very 
broad range is attractive for numerous applications, and in par- 
ticular in the field of transparent conductors. Carrier densi- 
ties > 10 20 cm" 3 , necessary for thin film operation, can be 
reached without extrinsic metal dopants, which are additional 
scattering centers [10]. Moreover, the possibility of after- 
growth patterning of conductive paths could yield important 
practical advantages. Namely, the initial growth conditions 
could be set independently of the required final conductiv- 
ity, and etching processes often involved in oxide structuring 
could potentially be avoided. 

Finally, the present study is likewise relevant for anatase- 
based devices employing nanostructured materials, where the 
overall transport properties depend on inter- as well as intra- 
particle processes. The crossover from a polaronic to a diffu- 
sive regime is expected to occur when the overlap between the 
polaron clouds becomes significant. The ARPES data of Fig. 
4 suggest that this happens around n* ~ 10 19 cm" 3 . Estimat- 
ing the polaron radius r p from the average separation between 
polarons d ~ n" 1 / 3 = 2r p , gives r p ~ 20 A. By comparison, 
Frohlich's model in the intermediate coupling regime yields 
r p = y/h/2mFujQ ~ 10 A. Both values are much smaller 
than the typical dimensions (few nm) of anatase nanoparticles 
considered for applications. Therefore, the polaronic nature of 
the QPs and their evolution upon electron doping must neces- 
sarily be taken into account when modeling transport in actual 
devices. 
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